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Kurzreferat 


Ziel der Arbeit ist es, die Diederich-Methode zur Berechnung der Auftriebsverteilung eines 
Tragflügels im Tabellenkalkulationsprogramm Microsoft Excel basierend auf didaktischen 
Überlegungen zur Verfügung zu stellen. Die Diederich-Methode wird basierend auf Primär- 
und Sekundärliteratur beschrieben. Diagramme werden digitalisiert, damit die Methode 
automatisch ablaufen kann. Zur Optimierung der Auftriebsverteilung des Flügels werden die 
elliptische und die dreieckige Auftriebsverteilung sowie die Auftriebsverteilung nach Mason 
zum Vergleich angeboten. Es wird eine Methode zur Berechnung des maximalen 
Auftriebsbeiwertes des Flügels in die Diederich-Methode integriert. Dazu müssen die 
maximalen Auftriebsbeiwerte der Profile an der Flügelwurzel und an der Flügelspitze in das 
Programm eingegeben werden. Die Berechnung setzt einen Trapezflügel voraus. Sowohl 
Flügelpfeilung als auch eine lineare Flügelverwindung können berücksichtigt werden. Die 
Streckung darf keine zu kleinen Werte annehmen. Es wird subsonische Strömung 
vorausgesetzt und eine Strömung ohne Ablösung. Da nur der Flügel beschrieben wird, bleiben 
alle weiteren Einflüsse wie beispielsweise vom Rumpf oder von den Triebwerken 
unberücksichtigt. Die Excel-Arbeitmappe wurde für die Lehre im Flugzeugvorentwurf 
erstellt. Derzeit wird die Diederich-Methode offenbar nirgends als Tabellenkalkulation 
angeboten. Mit dieser Arbeit kann diese Lücke geschlossen werden. 


Abstract 


Calculating the Wing Lift Distribution With the Diederich Method in Microsoft Excel 


Aim of this project is to provide the Diederich Method for calculating the lift distribution of a 
wing in a Microsoft Excel spreadsheet based on didactic considerations. The Diederich 
Method is described based on primary and secondary literature. Diagrams are digitized so that 
the method can run automatically. To optimize the lift distribution of the wing, the elliptical 
and triangular lift distribution as well as Mason's lift distribution are offered for comparison. 
A method for calculating the maximum lift coefficient of the wing is integrated into the 
Diederich Method. To do this, the maximum lift coefficients of the airfoils at the wing root 
and at the wing tip must be entered in the program. The calculation assumes a trapezoidal 
wing. Both wing sweep and linear wing twist can be taken into account. The aspect ratio must 
not assume values that are too small. Subsonic flow and unseparated flow are assumed. Since 
only the wing is described, all other influences such as from the fuselage or from the engines 
are not taken into account. The Excel workbook was created for teaching in aircraft 
preliminary design. At the moment, the Diederich Method is apparently nowhere offered as a 
spreadsheet. With this work, this gap can be closed. 
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Die Diederich-Methode zur Berechnung der 
Auftriebsverteilung am Tragflügel in 
Microsoft Excel 


Aufgabenstellung für ein Projekt 


Hintergrund 

Die Auftriebsverteilungen am Tragflügel ist von Bedeutung für den induzierten Widerstand, 
für Lasten am Flügel und somit für die Flügelmasse sowie für das Überziehverhalten und den 
maximalen Auftriebsbeiwert des Flügels. Bei der Diederich-Methode handelt es sich um eine 
semi-empirische Methode zur Bestimmung der Auftriebsverteilung von Tragflächen. Die 
Methode wurde beschrieben in DIEDERICH, Franklin W., 1952. A Simple Approximate 
Method for Calculating Spanwise Lift Distributions and Aerodynamic Influence Coefficients 
at Subsonic Speeds. Washington: NACA (Technical Note 2751). Einen leichteren Zugang 
findet man in TORENBEEK, Egbert, 1988. Synthesis of Subsonic Airplane Design. Delft: 
University Press. Verfügbar unter: https://bit.ly/3m8KIIV. Eine Excel-Arbeitsmappe mit der 
Methode wurde von Priyanka Barua und Dieter Scholz bereits erstellt, kann genutzt werden, 
bedarf aber noch einer didaktischen Überarbeitung. 


Aufgabe 

Die Diederich-Methode soll als Excel-Tabelle zur Verfügung gestellt werden, in einer Form 

angepaßt an vorhandene Werkzeuge zum Flugzeugentwurf an der HAW Hamburg. Dies sind 

die Arbeitsschritte: 

e Aufarbeitung der mathematischen Grundlagen der Methode. 

e Funktionale und optische Überarbeitung der vorhandenen Excel-Arbeitsmappe. 

e Erstellung einer Internetseite und einer Bedienungsanleitung zur Bereitstellung der 
Methode online in Englisch. 


Die Ergebnisse sollen in einem Bericht dokumentiert werden. Es sind die DIN-Normen zur 
Erstellung technisch-wissenschaftlicher Berichte zu beachten. 
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1 Einleitung 


1.1 Motivation und Ziel der Arbeit 


Der Entwurf eines neuen Flugzeugs umfasst eine Vielzahl von Themengebieten, die alle 
individuell einen direkten und wichtigen Einfluss, auf das am Ende fertig gestellte Flugzeug, 
ausüben. Eines dieser Entwurfsgebiete ist die Entwicklung der Tragfläche des Flugzeugs. Sie 
stellt ein wesentliches Element im Entwurfsablauf dar. Die Auslegung des Flügels 
entscheidet am Ende maßgebend über zentrale Eigenschaften wie Verbrauch, 
Flugeigenschaften und Gewicht des Flugzeugs. Um ein Vergleichsmaß für unterschiedliche 
Tragflächenentwürfe zur Verfügung zu haben, wird ein besonderes Augenmerk auf die 
Auftriebsverteilung gelegt. Sie gibt einen direkten Aufschluss über die Eigenschaften, die die 
Tragfläche aufweist. Unteranderem kann entnommen werden, wie die Strukturbelastung des 
Flügels selbst sowie die am Rumpfübergang erbrachte Biegemomentbeanspruchung ausfällt. 
Gleichermaßen lassen sich Aussagen über entstehende Widerstände am Tragflügel treffen, 
welche in Bezug auf die Wirtschaftlichkeit von großer Bedeutung sind. Da die Berechnung 
einer Auftriebsverteilung bei Berücksichtigung zahlreicher physikalischer Größen sowie 
Flügel und Atmosphärenparameter mit erhöhtem Aufwand verbunden ist, stellt sich die Frage, 
wie man diese Arbeit, bei Beibehaltung der Ergebnisgenauigkeit minimieren kann. Firmen 
verfügen an dieser Stelle über komplexe Rechenprogramme, die der Öffentlichkeit jedoch 
aufgrund von Geheimhaltung oder Kostengründen nicht zur Verfügung stehen. Bei Studenten 
besteht allerdings häufig ein großes Interesse die Vorgänge an einer Tragfläche zu verstehen 
und gegebenenfalls auf eigene Rechnungen oder Entwürfe anzuwenden. Dabei ist es 
wünschenswert alle Möglichkeiten kompakt und verständlich in einem Programm integriert 
zu haben. Ziel dieser Arbeit ist somit die Erstellung eines solchen Programms, welches unter 
Einhaltung eines didaktisch bedachten Aufbaus, eine möglichst intuitive Benutzung 
gewährleistet. Das Programm, welches in Form einer Excel-Arbeitsmappe zur Verfügung 
stehen wird, ist besonders an Studierende im Bereich Flugzeugbau gerichtet, steht aber 
selbstverständlich auch jedem anderen interessierten Nutzer zur Verfügung. 
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1.2  Begriffsdefinitionen 


Aerodynamik 
Die Aerodynamik beschreibt die Vorgänge, die bei der Umströmung von Körpern durch 
Gase, insbesondere durch Luft entstehen (Duden 2020). 


Diederich-Methode 
Die Diederich-Methode ist ein Berechnungsverfahren zur Beschreibung der 
Auftriebsverteilung über einer Tragfläche. 


Lehrmittel 
Lehrmittel sind Materialien einer Lehrveranstaltung, die unterstützend zur Erzielung eines 
Lernziels eingesetzt werden. 


1.3  Literaturübersicht 


Als wichtigste Literaturquelle dieser Arbeit ist Torenbeek 1988 aufzuführen. Torenbeek führt 
die Verfahren der Diederich-Methode, die durch Diederich 1952 publiziert wurden, erneut 
auf und kombiniert diese mit den Erkenntnissen aus Anderson 1936. Aufgrund der 
Zusammenführung und der einfachen Formelverwertbarkeit des nach Torenbeek 1988 
aufgeführten Ablaufes, bauen die mathematischen Grundlagen dieser Arbeit weitestgehend 
auf Torenbeek 1988 auf. Für die Oberflächengestaltung der Excel-Arbeitsmappe wurde 
Bartels 2003 herangezogen. Die Erstellung der Website beruht weitgehend auf den 
Grundlagen in Haunschild 2017. 


1.4 Aufbau der Arbeit 


Die Erstellung der Lehrmittel wird in den folgenden Abschnitten wie folgt bearbeitet: 


Abschnitt 2 befasst sich mit Auftriebsverteilungen im Allgemeinen und behandelt die 
Auswirkungen bestimmter Verteilungsformen. 


Abschnitt 3 behandelt die mathematischen Grundlagen zur Berechnung der 
Auftriebsverteilung nach der Diederich-Methode. 


Abschnitt 4 geht auf die Oberflächengestaltung eines Lehrmittels ein und beschreibt 
anschließend die Erstellung und Funktionsweise der Excel-Arbeitsmappe. 


Abschnitt 5 


Abschnitt 6 


Anhang A 


Anhang B 


Anhang C 
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erläutert den Aufbau und die Gestaltung der Webseite. 

fasst die erzielten Ergebnisse der Projektarbeit zusammen. 

beinhaltet die Gleichungen zur Bestimmung benötigter Atmosphärendaten. 
enthält die Quellcodes der in Visual Basic erstellten Funktionen. 


beinhaltet die Bedienungsanleitung der Excel-Arbeitsmappe. 
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2  Auftriebsverteilungen an Flugzeugtragflächen 


Unter einer Auftriebsverteilung versteht man die lokale Verteilung der am gesamten Flugzeug 
wirkenden Auftriebskräfte. Dabei wird allgemein der Auftrieb von Rumpf, Leitwerken sowie 
weiteren Flugzeugkomponenten berücksichtigt. Für die Vordimensionierung werden diese 
Faktoren häufig vernachlässigt und stattdessen eine Betrachtung der gesamten 
Grundrissfläche der Tragfläche, inklusiver der vom Rumpf bedecken Bereiche berücksichtigt 
(vgl. Torenbeek 1988). Die Auftriebsverteilung wird häufig mit dem griechischen 
Buchstaben Gamma 7’ bezeichnet. Sie kann über den zweidimensionalen Auftriebsbeiwert des 
Profils c, multipliziert mit der lokal vorliegenden Profiltiefe c(n) definiert werden. Um eine 
größenunabhängige Verteilung zu erhalten wird die Profiltiefe c(n) auf die mittlere Profiltiefe 
C, der Tragfläche bezogen. Die Flügelkoordinate 7 beschreibt die Tragfläche entlang der 
Spannweite. Die beinhalteten Größen werden in Abschnitt 3.1 noch einmal ausführlich 
behandelt. Die Auftriebsverteilung lässt sich nun beschreiben zu: 


= cıc 
I’n) = ne (m) (2.1) 
g 


2.1 Verteilungsformen und deren Bedeutung 


Auftriebsverteilungen können je nach Wahl der Flügelparameter unterschiedliche Formen 
annehmen. Die Form der Verteilung bestimmt dabei maßgebend die Auswirkungen auf die 
technischen Aspekte des Flugzeugs und sollte daher einer genauen Beobachtung unterliegen. 
Die nachfolgenden Verteilungsformen verfügen über gesonderte Eigenschaften, die hier kurz 
erläutert werden. 


2.1.1 Dreieckige Auftriebsverteilung 


Die dreieckige Auftriebsverteilung ist die einfachste Annahme einer Auftriebsfunktion. Sie ist 
häufig nicht der Realität entsprechend und findet sich in einer ähnlichen Form nur bei sehr 
spitz zulaufenden Flügeln mit hoher Streckung. Jedoch hat die dreieckige Verteilung, 
aufgrund der vergleichsweise rumpfnahen Ansammlung von Auftriebskräften, eine positive 
Auswirkung auf die Biegemomentbelastung an der Flügelwurzel. 
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2.1.2 Elliptische Auftriebsverteilung 


Die elliptische Auftriebsverteilung für alle Anstellwinkel wird bei einer elliptischen 
Profiltiefenverteilung sowie einer parallelen Anordnung und dem Vorliegen geometrisch 
gleicher Flügelschnitte erreicht (Dubs 1979). Ist der Verlauf elliptisch bedeutet dies, dass die 
örtliche Vertikalgeschwindigkeit an der Flügelhinterkante über die gesamte Spannweite einen 
konstanten Verlauf aufweist. Der aus den Abwinden der Tragfläche resultierende induzierte 
Widerstand erreicht für diesen Fall einen Minimalwert. Geringe induzierte Widerstände sind 
insbesondere aufgrund wirtschaftlicher Aspekte von großer Bedeutung, da der Verbrauch eins 
Flugzeugs von den erzeugten Widerständen abhängt. Die elliptische Verteilung gilt somit 
häufig, insbesondere für das Erreichen hoher Gleitzahlen, als die optimale 
Auftriebsverteilung. 


2.1.3 Auftriebsverteilung nach Mason 


Die Auftriebsverteilung nach Mason stellt eine Mischform der beiden zuvor erläuterten 
Verteilungen dar. Sie berücksichtigt die erhöhte Strukturbelastung einer Tragfläche mit 
elliptischer Auftriebsverteilung im Vergleich zu einer mit dreieckiger Auftriebsverteilung. Bei 
einer Anpassung der elliptischen Verteilung an die Dreiecksverteilung wird zwar der 
induzierte Widerstand erhöht, das Strukturgewicht hingegen sinkt, da weniger 
Biegemomentbelastung an der Flügelwurzel aufgenommen werden muss. Mason beschäftigte 
sich nun mit der Frage bis zu welchem Punkt die Gewichtsreduzierung der Struktur einen 
größeren Vorteil bietet als die gleichzeitige Erhöhung des induzierten Widerstandes Nachteile 
mit sich bringt. Die Ergebnisse sind in Mason 2001 zusammengetragen. Die 
Auftriebsverteilung nach Mason ist in Bild 2.1, mit einer Reduzierung der 
Biegemomentbelastung an der Flügelwurzel von 11 %, in Blau zu sehen. In Rot ist die 
elliptische Verteilung und in Grün die dreieckige Verteilung dargestellt. 
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Auftriebsverteilungen über der Halbspannweite 


















































je —— Triangle 
= 1,00 
0 ——Elliptic 
——- Mason 
0,50 
0,00 
{0} 01 02 013 04 015 06 07 08 019 1 
-0,50 
n=2y/b 
Bild 2.1 Auftriebsverteilung nach Mason mit dreieckiger und elliptischer Verteilung als 
Vergleich 


2.2 Möglichkeiten der Berechnung 


Für die Berechnung von Auftriebsverteilungen existieren mehrere Verfahren. Die wohl 
bekannteste und noch heute für den ersten Entwurf eines Flugzeugs angewandte Methode ist 
die Berechnung nach der Prandtlschen Tragflügeltheorie. Sie basiert auf der Annahme einer 
um den Tragflügel wirkenden Zirkulationsverteilung, die nach der Kutta Joukowski- 
Bedingung proportional zur Auftriebsverteilung der Tragfläche ist (Oertel 2017). Die 
Berechnung erfolgt nach der sogenannten Prandtlschen Integro-Differantialgleichung. Da sich 
die Lösung der Integro-Differentialgleichung mit einfachen Rechenmethoden als schwierig 
gestaltet, bietet sich als weiteres Verfahren zur Berechnung der Auftriebsverteilung die 
Diederich-Methode an. Sie basiert auf der Annahme die Auftriebsverteilung in zwei 
Komponenten aufzuteilen und diese gesondert zu berechnen. Am Ende werden beide 
Verteilungen zur gesamten Auftriebsverteilung superpositioniert. Aufgrund der besseren 
Integration in eine Excel-Arbeitsmappe und dem leichteren Verständnis des Verfahrens, was 
insbesondere für den Einsatz innerhalb eines Lehrmittels von Bedeutung ist, wird in den 
folgenden Abschnitten mit der Diederich-Methode weitergearbeitet. 
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3  Auftriebsverteilung nach Diederich 


3.1. Theoretische Grundlagen 


Um die Berechnungsverfahren der Diederich-Methode zu verstehen, ist es notwendig 
zunächst die wichtigsten theoretischen Grundlagen zu erläutern. Die Gleichungen stammen, 
sofern nicht anders gekennzeichnet, aus Torenbeek 1988. 


3.1.1 Flügelkoordinate 


Zur Beschreibung einer Halbspannweite wird in dieser Arbeit die auf die Spannweite b 
normierte Koordinate n verwendet und ist definiert zu: 


2.y 
= — 3.1 


Der Nullpunkt liegt damit auf der Rumpflängsachse. n beschreibt somit im Bereich von -1 bis 
1 die Spannweite einer Tragfläche. 


3.1.2 Flügelfläche 


Die Flügelfläche S ist die Grundrissfläche des Flügels und kann mit folgender Formel nach 
Scholz 2013 bestimmt werden. 


1 
sS=2 | c(n) dy (3.2) 
0 


Dabei ist c(n) die Profiltiefenverteilung entlang der Tragfläche. 
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3.1.3 Streckung 
Die Streckung A ist wie folgt definiert und beschreibt den Schlankheitsgrad einer Tragfläche. 


b? 
A= = (3.3) 


Eine große Streckung deutet auf einen langgezogenen schmalen Flügel hin, während eine 
kleine Streckung einen kurzen und breiten Flügel beschreibt. 


3.1.4 Zuspitzung 


Die Zuspitzung A beschreibt das Verhältnis aus Profiltiefe an der Flügelwurzel zu Profiltiefe 
an der Flügelspitze und ist demnach definiert zu: 


c 
= 3.4) 
Cr 


Kleine Zuspitzungen treten auf bei sehr spitz zulaufenden Tragflächen, während eine 


Zuspitzung A = 1 einem Rechteckflügel entspricht. Für A = 0,45 wird nach Raymer 1992 
eine fast optimal elliptische Verteilung des Auftriebs erreicht. 


3.1.5 Profiltiefenbezeichnungen 


Die mittlere geometrische Profiltiefe c, ist die Profiltiefe eines äquivalenten ungepfeilten 
Rechteckflügels, der die gleiche Grundfläche wie die Ausgangstragfläche besitzt. Sie lässt 
sich über folgenden Ansatz, hergeleitet aus Gleichung 3.3, ermitteln. 


Gl 3.5) 


Die Profiltiefe an der Flügelwurzel c, und and der Flügelspitze c, lassen sich bestimmen zu: 


"“ G.6) 





q=i:Cc (3.7) 
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Liegt eine lineare Verteilung der Profiltiefe über der Halbspannweite c(n) vor, so kann diese 
mit einer einfachen linearen Funktion wie folgt beschrieben werden: 


ce) le-)ntG (3.8) 


3.1.6 Flügelpfeilung 


Der Pfeilungswinkel @,; wird definiert zwischen der senkrechten zur Rumpflängsachse und 
der 25%-Linie. Die 25%-Linie verbindet die Punkte der einzelnen Profilschnitte, die je bei 
1/4 der Profiltiefe von der Profilnase ausgehend, liegen. Wird der Pfeilungswinkel der 
Tragfläche nicht mit der 25%-Linie angegeben so lässt er sich mit folgender Gleichung nach 
Scholz 2013 von einer m-Linie auf eine n-Linie umrechnen. Die Werte m und n werden 
jeweils in Prozent angegeben. 


2 , : 4 m-m 2) 3.9) 
Anm arcen anlem) Bauens | ' 


3.1.7. Verwindung 


Die Verwindung e beschreibt die Veränderung des Einstellwinkels über der Halbspannweite. 
Handelt es sich um eine negative Verwindung senkt sich die Profilnase zur Flügelspitze 
langsam ab. Hier spricht man von einem wash-out. Liegt eine positive Verdrehung vor, wird 
diese als wash-in bezeichnet. Bei der Beschreibung der Verwindung wird der Einstellwinkel 
der Wurzelrippe als Bezug gewählt. Dieser wird bei den folgenden Berechnungen mit null 
Grad vorausgesetzt. Somit folgt mithilfe beider Anstellwinkel &, und a, die Flügelspitzen- 
verwindung nach Torenbeek 1988 zu: 


& =. — 0, (3.10) 


Liegt eine lineare Verwindung über der Halbspannweite vor, so lässt sich diese in 
Abhängigkeit von n beschreiben zu: 


em)=&:'n (3.11) 
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3.1.8 Relative Profildicke 


Die relative Profildicke (t/c) ist das Verhältnis aus der dicksten Stelle eines Profils und der 
dazugehörigen Profiltiefe. Ist diese nicht vorhanden, so kann sie in Abbott 1959 für eine 
große Auswahl an Profilen eingesehen werden. Die relative Profildicke lässt sich hier an den 
letzten beiden Ziffern der Nennung erkennen und wird in Prozent angegeben. 


3.1.9 Machzahl 


Die dimensionslose Machzahl M beschreibt die Fluggeschwindigkeit bezogen auf die örtlich 
vorliegende Schallgeschwindigkeit. Die Schallgeschwindigkeit hängt dabei von der Flughöhe 
ab und lässt sich aus den Atmosphärengleichungen der ISA ermitteln, die in Anhang A nach 
Scholz 2012 aufgezeigt sind. Für die Machzahl folgt schließlich: 


V 
M= 


_ (3.12) 
a 


3.1.10 Auftriebsbeiwert 


Der Auftriebsbeiwert C, ist eine dimensionslose Zahl zur Beschreibung des 
Auftriebsverhalten einer Tragfläche. Ermittelt wird dieser nach Dubs 1979 über folgende 
Gleichung: 


2-L 


ne 
IT 0-V2.S 


(3.13) 


Dabei kann für den Horizontalflug der Auftrieb L dem Flugzeuggewicht gleichgestellt 
werden. V ist die Fluggeschwindigkeit wie auch bei der Machzahlermittlung im vorherigen 
Abschnitt und S die Flügelfläche aus Abschnitt 3.1.2. Die Dichte p kann ebenfalls über die 
Gleichungen in Anhang A ermittelt werden. 


3.1.11 Reynoldszahl 


Die Reynoldszahl Re ist eine dimensionslose Kenngröße der Aerodynamik, die unter anderem 
zur Unterscheidung von laminarer und turbulenter Strömung verwendet wird. Sie ist definiert 
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als das Verhältnis aus dem Produkt der Länge des umströmten Körpers ! und der 
Strömungsgeschwindigkeit V zur kinematischen Viskosität v. 


Ba (3.14) 
Vv 


Für den angewandten Fall an der Tragfläche kann die Länge des umströmten Körpers mit der 
mittleren Profiltiefe c, aus Abschnitt 3.1.5 gleichgesetzt werden. Die kinematische Viskosität 


v ist ebenfalls über die Gleichungen in Anhang A zu ermitteln. 


3.2 Berechnung der Auftriebsverteilung 


Bei der in dieser Arbeit verwendeten Methode zur Berechnung der Auftriebsverteilung 
handelt es sich um eine leicht modifizierte Variante, des in Diederich 1952 beschriebenen 
Verfahrens. Mit kleinen Einschüben wird im Grundsatz, der von Torenbeek 1988 aufgeführte 
Ablauf aufgrund der besseren Didaktik verwendet. Torenbeeks Arbeit basiert auf der von 
Diederich und kombiniert diese mit den Erkenntnissen von Anderson 1936. 


Bei der Diederich-Methode handelt es sich um eine semi-empirische Methode zur 
Bestimmung der Auftriebsverteilung von Tragflächen. Das Verfahren basiert auf der 
Annahme die Auftriebsverteilung in zwei Komponenten aufzuteilen. Dabei handelt es sich 
zum einen um die grundlegende Auftriebsverteilung, im Folgenden als „basic lift distribution“ 
bezeichnet und zum anderen um die zusätzliche Auftriebsverteilung, im Folgenden als 
„additional lift distribution“ bezeichnet. Die basic lift distribution ist definiert als 
Auftriebsverteilung für den verwundenen Flügel mit einem Auftriebsbeiwert von C, = 0. Die 
additional lift distribution ist die Auftriebsverteilung an dem unverwundenen Flügel mit dem 
normierten Auftriebsbeiwert von C, =1. Die Superposition beider Verteilungen und 
Anpassung an den vorhandenen Auftriebsbeiwert liefert schließlich die Auftriebsverteilung 
der Tragfläche. 


Damit die Methode zielführende Ergebnisse liefert müssen nach Torenbeek 1988 folgende 
Bedingungen erfüllt sein. 


e Die Fluggeschwindigkeit liegt innerhalb des subsonischen Bereiches. Nach Roskam 
1997 beginnt ab Machzahlen von M > 0,8 der transsonische Bereich. Um zu 
gewährleisten, dass keine oder vernachlässigbar kleine Kompressibilitätseffekte 
auftreten, sollte dieser Wert nicht überschritten werden. 

e Die Anstellwinkel sind verhältnismäßig klein, sodass die Strömung zu jeder Zeit 
anliegt. Ablösungsphänomene werden nicht berücksichtigt. 
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e Bei der Betrachtung der Tragfläche wird der Einfluss vom Rumpf unberücksichtigt 
gelassen. Es wird angenommen, dass die Eigenschaften des Rumpfes an relevanten 
Stellen, denen der alleinstehenden Tragfläche gleichgesetzt werden. 

e Die Verwindung der Tragfläche wird als linear vorausgesetzt. (siehe Abschnitt 3.1.7) 

e Die Zuspitzung des Flügels verläuft geradlinig. 

e Der Pfeilungswinkel liegt im Bereich —35° < 935 < 35°. 

e Die minimale Streckung A der Tragfläche lässt sich abhängig vom Pfeilungswinkel 
über die Gleichung 3.15 ermitteln. 


4 


A> — 3.15 
2 08 (9) So 


e Triebwerkseinflüsse, Bodeneffekte sowie Einwirkungen der Aeroelastizität bleiben 
unberücksichtigt. 


3.2.1 Additional Lift Distribution 


Für die Berechnung der additional lift distribution /,(n) wird zunächst die lift distribution 
function L. (n) nach Gleichung 3.16 berechnet. 


4 
um 4o.2. Mm +G-fo) 3.16) 
g 


Besteht der Sonderfall eines ungepfeilten Flügels (@35 = 0) äußert sich dies in einer 
elliptischen Verteilung der Pfeilungswinkelkorrekturfunktion f(n) siehe Bild 3.3. Folgend 
kann die additional lift distribution function vereinfacht werden zu: 


4 
= 44-2 Im 3.1) 


Cg 
Die benötigten Profiltiefen c(n) und c, sind bereits aus Abschnitt 3.1.5 bekannt. 


Die Koeffizienten C,, C2 und Cz lassen sich aus Bild 3.1 mithilfe des Grundrissparameters F 
bestimmen. Dieser ist definiert zu: 


Be 3.18) 
14 Fa c0os(P25) 


Der benötigte Auftriebsgradient c,, beschreibt den örtlichen Auftriebsanstieg einer 


Flügelsektion. Sollten keine genauen Daten zum Profil vorliegen kann nach Torenbeek 1988 
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der Wert c,, = 6,1 1/rad angenommen werden. Für relative Profildicken 0,1 < (t/c) < 


0,2 liefert dieser Wert zufriedenstellende Ergebnisse. 


Liegen Profildaten vor so kann c,,unter zusätzlicher Berücksichtigung des 
Machzahleinflusses nach DATCOM 1978 bestimmt werden. 


a le 3.19) 
re Ve ee 15 03 5 
ß theory 
Mit: 
t 1 
GHRTER = 27 +4,7:— (14 0,00375- P75) (3.20) 


Der Winkel der Profilhinterkante #7; kann wie in Bild 3.2 beschrieben ermittelt werden. Hier 
kann ebenfalls mithilfe der in Abschnitt 3.1.11 erläuterten Reynoldszahl und dem ermittelten 
Profilhinterkantenwinkel das benötigte Verhältnis aus dem experimentellen 


Auftriebsgradienten c,, zum theoretischen Auftriebsgradienten RE bestimmt 
[4 


werden. 
Der noch fehlende Prandtl-Glauert Faktor ß lässt sich mit der Machzahl M bestimmen. 


(3.21) 





09202 4 6:77 87107 EE120 7 14 


2rtA 
? C] COSA 2, 


Bild 3.1 Ch, RR und C, in Abhängigkeit des Grundrissparameters F 
(Torenbeek 1988) 
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Note: Interpolate by plotting 
2 


rt vs logıo Re 





Bild 3.2 Definition des Profilhinterkantenwinkels or; und Ermittlung 
des Auftriebsgradientenverhältnisses (DATCOM 1978) 


Die noch benötigte Pfeilungswinkelkorrekturfunktion f(n) ist in Bild 3.3 dargestellt. Die 
Wahl der Kurve hängt dabei vom effektiven Pfeilungswinkel @g ab, der sich über die 


Gleichung 3.22 bestimmen lässt. 


_ er) 
PB = arctan rn 


Da das Diagramm nur für bestimmte effektive Pfeilungswinkel eine Kurvendarstellung liefert, 


muss bei speziellen Winkeln eine Interpolation vorgenommen werden. 
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Bild 3.3 Pfeilungswinkelkorrekturfunktion f(n)(nach Diederich 1952 aus Torenbeek 1988) 
Agist in diese Arbeit 9; . 


Die lift distribution function L. (n) kann nun mit der Gleichung 3.16 oder 3.17 ermittelt 
werden. Da eine Normierung auf den Auftriebsbeiwert von C, = 1 vorliegt, muss die 
Funktion noch mit dem vorhandenen Auftriebsbeiwert der Tragfläche multipliziert werden. Es 
ergibt sich die additional lift distribution /,(n). 


Ca € 
In) = (2) mM)=LamC, (3.23) 


g 


3.2.2 Basic Lift Distribution 


Für die basic lift distribution 7,(n) wird zunächst die basic lift distribution function L, (n) 
ermittelt. 


L, m) =La(m)'Cy:cos(Ag)- = + “) -ß-E (3.24) 
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Der Anteil e(n)/e; reduziert sich dabei aufgrund der Annahme einer linearen Verwindung mit 
der Flügelwurzel als Bezugsnullpunkt zu 





Em er. n (3.25) 
& & 


Der Koeffizient C, kann ebenfalls aus Bild 3.1 mithilfe des nach Gleichung 3.18 bestimmten 
Grundrissparameters F ermittelt werden. 


Der Jones Edges Velocity Faktor E kann mit folgender Näherung beschrieben werden zu: 
E=1+— — (3.26) 
Mithilfe der additional lift distribution function nach Gleichung 3.16 oder 3.17 und der 


Annahme einer linearen Verwindung wird der Faktor @,, zu 


1 
0, =- | n-La man (327) 
0 


Für ungepfeilte Flügel lässt sich die Gleichung ähnlich der Annahme für die Gleichung 3.17 


vereinfachen zu 


1:24 4 
le) (3.28) 


=D 3ER) 37 


=(] 


1 


Die basic lift distribution /,(n) lässt sich nun aus der basic lift distribution function L, (n) 
und den bereits ermittelten Faktoren berechnen. 


e,„.e L,m-&: u 
Tr u ee er ur 3.29 
pn) | , )m E (3.29) 


3.2.3 Gesamte Auftriebsverteilung 


Um nun die gesamte Auftriebsverteilung darzustellen werden die additional- und die basic lift 
distribution aus Gleichung 3.23 und Gleichung 3.29 addiert. 


TD)=-Rm)+NnG@m (3.30) 
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3.3 Verteilung der lokalen Auftriebsbeiwerte 


Die Verteilung der Auftriebsbeiwerte c,(n) ist unteranderem wichtig, um zu lokalisieren wo 
entlang der Tragfläche in kritischen Fluglagen die Strömung als erstes ablöst. Um diese zu 
ermitteln liegt nach Torenbeek 1988 folgender Zusammenhang vor. 


rn m \m+(fe 


Daraus lässt sich mithilfe der Gleichungen 3.23 und 3.29 folgender Zusammenhang 








en ) (3.31) 


formulieren. 


cam) = (Em - Lam: Cı 39) 


I»m-&: Ey 


r (3.33) 


c 

m=()m- 
Aus den bestimmten Größen kann nun die Verteilung der lokalen Auftriebsbeiwerte ermittelt 
werden. 


cm) = can) + c, Mm) (3.34) 


3.4 Maximaler Auftriebsbeiwert 


Die Berechnung des maximalen Auftriebsbeiwertes der Tragfläche erfolgt in Anlehnung an 
das in Torenbeek 1988 beschriebene Verfahren. Dieser wird erreicht, wenn an einer Stelle 
entlang der Tragfläche der lokale Auftriebsbeiwert des Profils den zu dieser Sektion 
zugehörigen maximalen Auftriebsbeiwert übersteigt. In Bild 3.4 ist dieser Moment 
beispielhaft dargestellt. Bei weiterer Erhöhung des Anstellwinkels beginnt die Strömung an 
dieser Stelle abzulösen und die Auftriebsleistung der Tragfläche sinkt. Mit den in den 
vorherigen Abschnitten ermittelten Größen kann der maximale Auftriebsbeiwert der 
Tragfläche C, max aus dem Minimum des eingeklammerten Terms und dem vorhandenen 
Auftriebsbeiwert der Tragfläche C, bestimmt werden. 


CLmax >= min 


een) LE I. ©; (3.35) 
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Der maximale Auftriebsbeiwert der Tragfläche wird erreicht, wenn an (mindestens) 
einer spannweitigen Position am Flügel der lokale maximale Auftriebsbeiwert erreicht 
wird. 


29 


4 Erstellung der Excel-Arbeitsmappe 


4.1 Allgemein 


Eine Möglichkeit, um Studierenden die Inhalte des Tragflächenentwurfs näher zu bringen, 
ohne sich intensiv mit der Literatur auseinanderzusetzen zu müssen, stellt die Bereitstellung 
einer Excel-Arbeitsmappe dar. Hier können über Parameteränderungen schnell die 
grundlegenden Prinzipien verstanden werden und die Thematik kann „spielerisch“ erlernt 
werden. Diese Arbeit baut auf einer bereits bestehenden Excel-Arbeitsmappe auf, die aus 
einer Zusammenarbeit von Prof. Dr.-Ing. Dieter Scholz und Priyanka Barua entstanden ist. 
Die grundlegenden Rechenvorgänge für die Auftriebsverteilung sind in dieser Datei bereits 
vorhanden gewesen. Damit das Lehrmittel für Studenten ohne direkte Fachkenntnisse nutzbar 
wird, wurden einige Ergänzungen sowie eine grundlegende didaktische Aufarbeitung 
vorgenommen. Außerdem wurde der gesamte Umfang der Lehrmittel in englischer Sprache 
erstellt, um auch nicht deutschsprachigen Studenten die Benutzung zu ermöglichen. Die 
nachfolgenden Abschnitte erklären die Funktionsweise sowie die Grundlagen zur 
Oberflächengestaltung der Excel-Arbeitsmappe. 


4.2  Gestaltungsgrundlagen 


Ein wichtiger Faktor, der zur Benutzerfreundlichkeit eines Lerntools beiträgt, ist die 
zielgerichtete Gestaltung der Benutzeroberfläche. Es ist wichtig, dass auf den ersten Blick 
keine Überforderung des Benutzers stattfindet, sondern dass intuitiv und ohne viel Einlesezeit 
klar ist, wie das Lehrmittel zu bedienen ist. Um diesen Effekt zu erreichen, ist es wichtig 
einige Grundregeln zu befolgen, auch wenn am Ende der eigene Eindruck ein wichtiges 
Kriterium darstellt. Ein bedeutender Faktor ist die Farbgebung der Oberfläche, die im 
Folgenden genauer erläutert werden soll. 


Die Farbgebung ist sehr vielfältig und richtet sich unteranderem nach der Zielgruppe, die das 
Lehrmittel voraussichtlich am häufigsten anwendet. Da hier keine Reize, sondern Struktur 
und Übersichtlichkeit die bestimmenden Faktoren sind, sollten beruhigende Farben gewählt 
werden. Hierfür eignen sich die Farben Blau und Grün. Sie stehen nach Bartel 2003 in 
kombinierter Anwendung für wenig Aufregung und Ruhe und sollen noch um die Farbe Weiß 
ergänzt werden. Da Weiß auf Computerbildschirmen häufig einen grellen Eindruck erzeugt, 
wurde das Weiß durch die leicht abgedunkelte Variante des Hellgraus ersetzt. Wichtig ist 
ebenso, dass sich die Farbwahl auf wenige Farben beschränkt. Buntheit senkt nicht nur die 
Seriosität, sondern gilt nach Bartel 2003 als physikalischer Reiz und sollte für ein Lehrmittel 
vermieden werden. Die gewählten Farben sollten sich außerdem über die gesamte 
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Benutzeroberfläche erstrecken, um ein einheitliches Gesamtbild zu erzeugen. Bei 
Kurvendarstellungen innerhalb eines Diagramms kann auf eine bunte Farbgebung jedoch 
nicht verzichtet werden, um die Differenzierbarkeit der Funktionen zu gewährleisten. Um 
bestimmte Aspekte hervorzuheben, eignet sich die Farbe Rot. Sie hat nach Bartel 2003 keine 
Tiefenwirkung und drängt sich somit in den Vordergrund des Betrachters. Bei der 
Anwendung sollte hierbei jedoch auf eine minimalistische Auswahl der Rotelemente geachtet 
werden, um den Benutzer nicht mit zu vielen Reizpunkten zu überfordern. 


Bei der Benutzeroberfläche der Excel-Arbeitsmappe wird folglich auf die Farben Blau und 
Grün als Hauptelemente gesetzt, um die Oberfläche mit einer Struktur zu versehen. Als 
Hintergrund dient das gewählte Hellgrau, welches farblos aber nicht ganz so aufdringlich wie 
Weiß wirkt. Die Eingabedaten sind rot gefärbt und vermitteln dem Benutzer sofort, wo 
Änderungen vorgenommen werden können und wo nicht. Untermauert wird dies durch die 
blaue Farbe der Eingabefelder und die grüne Farbe der Ausgabefelder. Der Benutzer 
entwickelt hierdurch ein direktes Verständnis wo nach Ergebnissen gesucht werden muss und 
wo die veränderlichen Parameter zu finden sind. Die rote Farbe lenkt den Blick schließlich 
nochmal auf das innerhalb des Eingabebereiches wichtige Zahlenfeld. 


Bei der Gestaltung der Oberfläche ist ebenfalls auf eine Gruppierung der Einzelelemente zu 
achten. Strukturen, die sich in Gruppen befinden können auf den ersten Blick 
zusammengefasst werden und reduzieren damit die für den Nutzer zu verarbeitenden 
Informationen. Die Farbgebung kann hier unterstützend eingesetzt werden. 
Zusammenfassende Überschriften ermöglichen dem Benutzer schließlich eine gedankliche 
Zuordnung der von den Gruppen übernommenen Aufgaben und gewährleisten so eine 
intuitive Benutzung. 


Die Anwendung der Gestaltungsregeln ist in Bild 4.1 in Form der Benutzeroberfläche der 
Excel-Arbeitsmappe zu sehen. 
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Bild 4.1 Benutzeroberfläche der Excel-Arbeitsmappe 
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4.3 Benutzung 


Folgend wird auf die Benutzung sowie die Funktionsweise der Excel-Arbeitsmappe 
eingegangen. Da die theoretischen Grundlagen bereits in Abschnitt 3 behandelt wurden, wird 
bei den Berechnungsbezeichnungen auf dieses verwiesen. Die Excel-Arbeitsmappe besteht 
aus mehreren Arbeitsblättern. Die für den Benutzer hauptsächlich interessante Seite ist das 
„Diederich‘s Method“ Arbeitsblatt. Die weiteren Arbeitsblätter erläutern die Funktionsweise 
der Hauptseite oder liefern Zwischenergebnisse, die dadurch gesondert eingesehen werden 
können. Diese sind für die Benutzung jedoch nicht erforderlich und stellen nur eine weitere 
Informationsquelle dar. Zusätzliche Informationen zum Verständnis der Excel-Arbeitsmappe 
sind weitläufig an Feldern mit roten Dreiecken durch Notizen hinterlegt. 


4.3.1 Diederich’s Method 


Auf dem „Diederich’s Method“ Arbeitsblatt (zu sehen in Bild 4.1) befinden sich die 
Eingabefelder in Blau sowie die Ausgabefelder in Grün auf der linken Seite. Das mittig 
liegende Diagramm ist in zwei Bereiche aufgeteilt. Der obere Bereich zeigt die 
Auftriebsverteilungen über einer Halbspannweite. Die jeweilige Kurve ist rechts neben dem 
Diagramm in der Legende namentlich zugewiesen. Unter den Auftriebsverteilungen befindet 
sich eine Visualisierung der den Parametern entsprechenden Halbspannweite. Auf der rechten 
Seite sind noch einmal kurz die wichtigsten Fakten zur Benutzung erläutert. So kann auch 
ohne ausführliche Anleitung, die später in Form einer Datei im PDF zur Verfügung gestellt 
wird, eine korrekte Benutzung sichergestellt werden. Der Inhalt dieses Bereiches wird hier 
nicht näher erläutert, da es sich um eine Zusammenfassung der wichtigsten Inhalte der 
folgenden Abschnitte handelt. Weiter unten auf der Seite findet sich eine große Wertetabelle, 
die die Plotwerte aller Graphen und Hilfsgrößen beinhaltet. Sie ist außerhalb des Sichtfeldes 
positioniert, da sie für die unmittelbare Benutzung nicht notwendig ist. Sollten jedoch 
bestimmt Zahlenwerte benötigt oder ein allgemein besseres Verständnis der Rechnung 
gewünscht werden, kann diese jederzeit eingesehen werden. 


Input 


Die Eingabefelder sind aufgeteilt in General und Additional Input. Die Sparte General Input 
beinhaltet alle für die Ermittlung der Auftriebsverteilung benötigten Werte. Der Additional 
Input liefert zusätzliche Anpassungen, die die Verteilungen noch besser auf die 
Flugzeugkonfiguration anpassen aber nicht unbedingt notwendig sind. Die hier beim Öffnen 
vorliegenden Werte sind stellvertretend für ein beispielhaft gewähltes Standardprofil der 
NACA-6-Serie, dem NACA 63,-615. Im Ausgangszustand wird außerdem eine Flughöhe h 
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von 10 km und eine mittlere Profiltiefe c, von | m angenommen. Hiermit werden qualitativ 
aussagekräftige Werte der Auftriebsverteilungen erreicht. Es wird beim Vorliegen 
detaillierterer Daten aber unbedingt zum Anpassen des Additional Inputs auf den spezifischen 
Fall geraten. 


Die Eingabeparameter des General Inputs richten sich im Grundsatz nach den in Abschnitt 
3.2 erwähnten Richtlinien. Die genannten Vorgaben für den Pfeilungswinkel 35 und die 
Machzahl M sind rechts neben den Diagrammen im „How to use“ Feld übersichtlich 
aufgelistet. Die dort aufgezeigte minimale Streckung A ändert sich mithilfe der Gleichung 
3.15 je nach Eingabe des Pfeilungswinkels @,;. Die Zuspitzung A der Tragfläche sollte für 
konventionelle Flügelkonfigurationen im Bereich 0 <A < 1 liegen. Die ersten drei Parameter 
des General Inputs sind hauptsächlich für das Aussehen des Flügelgrundrisses verantwortlich. 
Wie einzelne Änderungen dieser sich auf die Tragfläche auswirken, kann somit sehr gut im 
unteren Bereich des Diagramms nachgesehen werden. Die Winkeleingabe der Verwindung an 
der Flügelspitze g; sollte realistische Werte umfassen, die eine Strömungsablösung am Profil 
ausschließen. Der Auftriebsbeiwert der Tragfläche C, sollte bereits bekannt sein. Eine 
Berechnung erfolgt anderenfalls mithilfe der Gleichung 3.13. Eine Darstellung des General 
Inputs ist in Bild 4.2 zu sehen. 





General Input 














Aspect Ratio A 7 e 
Taper Ratio A 0,3 - 
Quarter Chord Sweep Po.25 25 deg 
Mach No M 0,7 - 
Aerodynamic twist atthe tip & -3 deg 
Lift coefficient C ıl - 
Bild 4.2 Eingabeparameter des General Inputs 


Der erste Parameter des Additional Inputs ist die Flughöhe h. Eingabewerte bis 20 km 
liefern korrekte Atmosphärendaten, die für die Bestimmung der Reynoldszahl nach Gleichung 
3.14 erforderlich sind. Der zweite Parameter der mittleren Profiltiefe c, hat ebenfalls 
Einwirkung auf die Reynoldszahl der Tragfläche. Die relativen Profildicken (t/c) sowie die 
Hinterkantenwinkel @7; für Flügelspitze und Flügelwurzel können dem NACA-Katalog in 
Abbott 1959 entnommen werden. Bei unterschiedlicher Verteilung über der Halbspannweite 
werden die Werte dementsprechend angepasst. Bei einem konstanten Profil des Flügels 
werden jeweils die gleichen Werte für Flügelspitze und Flügelwurzel eingegeben. In 
Abhängigkeit dieser Werte wird schließlich, wie in Abschnitt 3.2.1 behandelt, der 
Auftriebsgradient c,, einer Flügelsektion bestimmt. Der maximale Auftriebsbeiwert des 
Profils kann ebenfalls dem NACA-Katalog in Abbott 1959 entnommen werden. Wie bereits 
bei den vorherigen Werten können hier entweder gleiche oder bei verschiedenen Profilen 
unterschiedliche Werte für Flügelspitze und Flügelwurzel eingegeben werden. Nach dem 
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Verfahren in Abschnitt 3.3 wird hieraus der maximale Auftriebsbeiwert C, max der Tragfläche 
bestimmt. Die Parameter des Additional Inputs sind in Bild 4.3 dargestellt. 





Additional Input 











Mean geometric chord G; il m 
Flight altitude h 10 km 
Thickness ratio atthe root (t/e)- 15% - 
Thickness ratio atthe tip (t/e): 15% - 
Airfoil trailing edge angle atthe root Örer 14,00 deg 
Airfoil trailing edge angle at the tip Örgı 14,00 deg 
Maximum section lift atthe root sur 1,6 - 
Maximum section lift at the tip Cm 1,6 - 
Bild 4.3 Eingabeparameter des Additional Inputs 

Output 


Die Ausgabefelder sind aufgeteilt in drei Bereiche. Den „Comparison Output“, den 
„Additional Output“ sowie den „Calculated Values‘“ Output. 


Der Comparison Output liefert ein direktes Vergleichsmaß zwischen der nach der 
Diederich-Methode ermittelten Auftriebsverteilung und den in Abschnitt 2.1 behandelten 
Vergleichsverteilungen. Die Abweichung wird über die Summe der Residuenquadrate nach 
Barot 2020 innerhalb der Wertetabelle ermittelt. Ein kleiner Wert in der Ausgabe steht dabei 
für einen gut angenäherten Kurvenverlauf. Große Werte stehen hingegen für große 
Abweichungen. Der kleinste der drei Werte bestimmt somit die naheliegendste Eigenschaft, 
die die ermittelte Auftriebsverteilung aufweist. Die Eigenschaften der jeweiligen Verteilungen 
können in Abschnitt 2.1 eingesehen werden. In Bild 4.4 ist der Comparison Output 
dargestellt. 


Diederichs Lift Distribution Compared 


Error to Mason's lift distribution ee) 
| Er rorto tr ER RT + distribu ir n 





ngula 





Bild 4.4 Comparison Output mit Abweichungsmaß 


Der Additional Output liefert die für den Benutzer potenziell interessanten Zwischen- und 
Zusatzergebnisse. Unteranderem werden in Abhängigkeit der mittleren Profiltiefe c, die 
Profiltiefen der Flügelspitze c, und der Flügelwurzel c; für die eingegebenen Parameter des 
General Inputs ausgegeben. Weiterhin kann der maximale Auftriebsbeiwert C, max der 
Tragfläche und der Mittelwert des Auftriebsgradienten c,, aus den Auftriebsgradienten der 
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Flügelspitze c,,, und der Flügelwurzel c,,, abgelesen werden. Für die Ermittlung der 
Profildaten aus den NACA-Katalogen ist die Reynoldszahl der Tragfläche notwendig. Diese 
wird in Abhängigkeit von der Flughöhe h und der mittleren Profiltiefe c, ausgegeben und 
kann somit für die Profildaten im Additional Input verwendet werden. Für einige 
Berechnungen, insbesondere im Bereich des Flugzeugentwurfs, werden die Winkel der 
Vorderkante und der Hinterkante der Tragfläche benötigt. Diese können für den gewählten 
Pfeilungswinkel der 25%-Linie ebenfalls abgelesen werden. Weiterhin werden die 
Spannweite b und die Flügelfläche S sowie die Fluggeschwindigkeit in m/s und km/h 
angegeben. Die Auflistung des Additional Outputs ist in Bild 4.5 zu sehen. 


[Xefelhate]s 1 Keltigeltis 


RER ikas: 


m 


Maximum lift coefficient of the wing 
| t | 


Reynolds number Re 5.598.635,19 - 
Trailing edge angle i 13,25 deg 
Wing surface area 2.005 =) m? 


True airspeed Vras 754,63 km/h 
Bild 4.5 Ausgabewerte des Additional Outputs 


‚4 
‚40 B 
7,00 

















Der Calculated Values Output stellt zahlreiche Zwischenergebnisse der Berechnung der 
Auftriebsverteilungen gebündelt dar. Hier können Ergebnisse nachverfolgt werden und für 
den Fall einer händischen Berechnung stehen diese für Zwischenwerte zur Verfügung. Die 
Ergebnisse finden sich zum Großteil auch in den zusätzlichen Arbeitsblättern mit den 
jeweiligen Erzeugnisgraphen wieder. Hier kann noch einmal eingesehen werden, von wo die 
Ergebnisse stammen. Der Calculated Values Output ist in Bild 4.6 dargestellt. 


36 


[er [at] E12 RVETT-R 


2nA/ Cu C0Sg0.25 


Diederich's Factor 


Diederich's Factor 


Factor 








Effective sweep angle 


Lift curve slope of the airfoil atthe root Gier 8,37 1/rad 
Lift curve slope of the airfoil atthe tip Cat 8,37 1/rad 





Speed of sound in flight altitude a 299,46 m/s 
Bild 4.6 Ausgabewerte des Calculated Values Outputs 





Auftriebsverteilungsdiagramm über der Halbspannweite 


Das Auftriebsverteilungsdiagramm beinhaltet zahlreiche Kurvenverläufe, die über die in 
Abschnitt 3.1.1 definierte Flügelkoordinate n aufgetragen sind und einen direkten Vergleich 
untereinander ermöglichen. Dabei sind die rötlich gefärbten Kurven die Auftriebsverteilung 
nach der Diederich-Methode. Die gesamte Auftriebsverteilung, im Diagramm mit Gamma 
bezeichnet, ergibt sich wie bereits in Abschnitt 3.2 behandelt aus der additional lift 
distribution Gamma_a und der basic lift distribution Gamma_b. Die lokalen 
Auftriebsbeiwerte entlang der Halbspannweite c,(m), cıa(m) und cı,(n) sind durch die 
blautonigen Kurven dargestellt. Sie sind wie in Gleichung 3.31 ersichtlich unter anderem die 
Eingangsgrößen für die Auftriebsverteilungen. Hieraus ist erkenntlich, dass die 
Auftriebsverteilung die Profiltiefe über der Halbspannweite berücksichtigt. Bei einer 
Zuspitzung von A = 1 und damit einer konstanten Verteilung der Profiltiefen sind die Kurven 
T(n)und c,(n) identisch. Weiterhin im Diagramm zu erkennen ist die Verteilung des 
maximalen lokalen Auftriebsbeiwertes c] max über der Halbspannweite. Die Kurve ist in grün 
dargestellt und durch eine lineare Funktion, abhängig von den Parametern im Additional 
Input, für die maximalen Auftriebsbeiwerte an Flügelwurzel und Flügelspitze, bestimmt. Die 
grautonigen Kurven repräsentieren die Vergleichs-auftriebsverteilungen, die im Comparison 
Output gelistet sind. Hier kann noch einmal direkt überprüft werden an welchen Stellen die 
Abweichungen auftreten und wie man diese durch mögliche Konfiguration der Tragfläche an 
die gewünschte Verteilung anpassen kann. Mit einem mittigen Klick auf das Diagramm 
erscheinen rechts neben diesem drei Optionen. Wählt man das untere Feld an, können 
einzelne Kurven aus- und abgewählt werden, die im Diagramm sichtbar gemacht werden 
sollen. Dies bietet die Möglichkeit eines direkten Vergleichs zweier Kurven untereinander. 
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Die Halbspannweite unter dem Auftriebsverteilungsdiagramm ist wie bereits erwähnt eine 
Visualisierung der ersten drei General Input Parameter. Die dunkelblauen Linien bilden dabei 
den Grundriss der halben Tragfläche. Die hellblaue Linie stellt die 25%-Linie des Flügels dar 
an welcher auch der Pfeilungswinkel gemessen wird. Für ein leichteres Verständnis und ein 
direkteren Sachbezug kann hier im Vergleich mit den Auftriebsverteilungen herausgefunden 
werden welche Einwirkungen die Flügelparameter auf die Kurvenverläufe haben. Der 
Überblick über die Diagramme kann erneut in Bild 4.1 eingesehen werden. 


4.3.2 Mason's Curve 


Das „Mason's Curve“ Arbeitsblatt dient im Wesentlichen zur Erstellung der Kurve nach 
Mason im Auftriebsverteilungsdiagramm auf der Hauptseite. Hier wurde die bereits in 
Abschnitt 2.1.3 behandelte Auftriebsverteilung nach Mason 2001 über eine Wertetabelle 
geplottet und mit einer anschließenden Annäherung durch eine Polynomfunktion dritten 
Grades für Excel verwertbar gemacht. Die Wertetabelle befindet sich, wie auch bei den 
anderen Arbeitsblättern, im unteren Bereich der Seite. Um den direkten Vergleich zwischen 
den Auftriebsverteilungen auf der Hauptseite gewährleiten zu können, musste noch der 
Maßstab zwischen der Verteilung nach Diederich und der Verteilung nach Mason definiert 
werden. Als Anhaltspunkt hierfür dient die optimal elliptische Verteilung im Diagramm nach 
Mason, wie in Bild 4.7 zu sehen ist, die der Pfeilungswinkelkorrekturfunktion f(n) für einen 
effektiven Pfeilungswinkel @g = 0 gleichgestellt werden kann. Der Quotient der 
Achsenabschnitte beider elliptischer Verteilungen definiert den Maßstabsfaktor. Dieser ist im 
grün gefärbten Ausgabefeld dargestellt. Zusätzlich zum Maßstabsfaktor wird das Mason- 
Polynom noch mit dem Auftriebsbeiwert C, der Tragfläche multipliziert, da wie auch bei den 
anderen Verteilungsfunktionen eine Normierung auf C, = 1 vorliegt. 
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Bild 4.7 Arbeitsblatt „Mason's Curve“ (mit Diagramm aus Mason 2001) 
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4.3.3 Sweep Correction Function 


Das Arbeitsblatt ,„Sweep Correction Function“ ist eine Rekonstruktion der 
Pfeilungswinkelkorrekturfunktion f(n) nach Diederich 1952, die bereits in Abschnitt 3.2.1 
in Bild 3.3 dargestellt ist. Diese wurde der Genauigkeit halber über Polynome 6. Grades 
angenähert. Hier wird abhängig vom Eingangswert des effektiven Pfeilungswinkels og nach 
Gl. 3.22 die Funktion f(n) an die Wertetabelle der Hauptseite weitergegeben. Bei effektiven 
Pfeilungswinkeln, die zwischen den gegebenen Funktionen liegen, interpoliert, die in Visual 
Basic erstellte Funktion, die benötigten Ausgabewerte. Der Quellcode zur Funktion kann im 
Anhang B eingesehen werden. Die Funktion für den Pfeilungswinkel 9, = 0 stellt außerdem 
die Grundlage der optimal elliptischen Verteilung im Auftriebsverteilungsdiagramm der 
Hauptseite dar und wird, wie die Mason-Kurve, zur Anpassung mit dem Auftriebsbeiwert 
C, der Tragfläche multipliziert. Die Oberfläche des Arbeitsblattes ist in Bild 4.8 dargestellt. 
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Bild 4.8 Arbeitsblatt „Sweep Correction Function“ 


4.3.4 Factors of Diederich's Method 


Das „Factors of Diederich‘s Method“ Arbeitsblatt, dargestellt in Bild 4.9, dient zur Ausgabe 
der Diederich Faktoren C, bis C,, die in Abschnitt 3.2 behandelt wurden. Über Polynome 3. 
Grades wurden die Kurven aus Bild 3.1 näherungsweise abgebildet und in Visual Basic 
übertragen. Sie geben in Abhängigkeit des Grundrissparameters F die benötigten Faktoren 
aus. Die Ergebnisse für die gewählten Input-Parameter sind hier aufgeführt, finden sich aber 
auch im Calculated Values Output auf der Hauptseite wieder. Der Quellcode der Funktion 
kann in Anhang B eingesehen werden. 
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LLICHUERTEICHIEITER 
LE 5,798936223 
c1 | 0,292977827 
c2 0,455246672 








c3 | 0,24618788| 
a2 C4 | 0,407975104 
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Poly. (C3) 
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Bild 4.9 Arbeitsblatt „Factors of Diederich's Method“ 


4.3.5 Lift Curve Slope 


In diesem Arbeitsblatt wird abhängig von dem auf der Hauptseite gewählten 
Profilhinterkantenwinkel 97; ein dazugehöriges Verhältnis aus experimentellem Auftriebs- 
gradienten C,„; Zu theoretischem Auftriebsgradienten (Cj4)ineory ausgegeben. Dieses Verhältnis 
wird für die Berechnung des Auftriebsgradienten c,„ des Profils benötigt, dessen Berechnung 
bereits in Abschnitt 3.2.1 durch die Gleichungen 3.19 und 3.20 aufgezeigt wurde. Die 
Ausgabe erfolgt wie in Bild 3.2 dargestellt in Abhängigkeit von der Reynoldszahl. Die 
Näherungspolynome der drei Kurvenverläufe für die Reynoldszahlen 10°, 10’und 
10® wurden in Visual Basic mit einer Interpolationsfunktion verknüpft und liefern auch für 
Reynoldszahlen zwischen den gegebenen Kurvenverläufen aussagekräftige Werte. Die 
Berechnung erfolgt hier jeweils für den Auftriebsgradienten der Flügelspitze c,„, sowie für 
den Auftriebsgradienten der Flügelwurzel cj,,. Für die Weiterrechnung innerhalb der 
Diederich-Methode wird der Mittelwert der Auftriebsgradienten des Profils gebildet und im 
Additional Output ausgegeben. Die Oberfläche des Arbeitsblattes ist in Bild 4.10 dargestellt. 
Der Quellcode für diese Funktion kann in Anhang B eingesehen werden. 
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Effect of Trailing Edge Angle on Section Lift Curve Slope 
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Bild 4.10 Arbeitsblatt „Lift Curve Slope“ 


4.3.6 Atmosphere Parameters 


Das Arbeitsblatt der „Atmosphere Parameters“ stellt die Grundlage für die Berechnung der 
Reynoldszahl in Abhängigkeit von der Flughöhe dar. Für Flughöhen h zwischen O0 und 20 km 
werden die erforderlichen Daten ausgegeben und visuell durch Diagramme unterstützt. Die 
Berechnung und die Erstellung der Diagramme, zu sehen in Bild 4.11, wurden mit den 
Gleichungen in Anhang A durchgeführt. Die Zwischenergebnisse für die auf der Hauptseite 
ausgewählte Höhe werden rechts im grünen Ausgabefeld aufgezeigt, befinden sich aber auch 
im Calculated Values Output auf der Hauptseite. 


Kinematic Viscosity Speed of Sound 























Bild 4.11 Arbeitsblatt „Atmosphere Parameters“ 
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5  Webseitenerstellung 


Um Interessierten Benutzern die Excel-Arbeitsmappe zur Verfügung stellen zu können, eignet 
sich die Erstellung einer Webseite. Hier kann per Mausklick ein sofortiger Zugang zur Datei 
sowie zur Bedienungsanleitung gewährleistet werden. Gleichzeitig bietet sie einen Ort an, der 
immer wieder zugänglich ist, um die gewünschte Datei zu erreichen. So kann auch bei Verlust 
von Daten oder dem fehlerhaften Bearbeiten der Excel-Arbeitsmappe eine neue Kopie auf den 
Computer geladen werden. 


5.1 Grundidee der Gestaltung 


Die Gestaltung der Webseite sollte nach einem einfachen Aufbau erfolgen. Dies liegt 
unteranderem daran, dass diese nur die Aufgabe der Datenüberbringung besitzt. Sie soll in 
keiner Form einen Nutzer bewerben oder von etwas überzeugen. Ein überladene Internetseite 
kann den Nutzer beim Aufruf durch Reizüberflutung davon abbringen sich mit der Thematik 
zu beschäftigen. Die Webseite bündelt somit die Interessenpunkte auf eine kompakte 
Darstellung. Diese beinhalten selbstverständlich die Überschrift, um einen Sachbezug 
herzustellen sowie ein Bild, das den sofortigen Einblick in die Benutzeroberfläche 
gewährleistet. Weiterhin werden die Excel-Arbeitsmappe und die Bedienungsanleitung 
gesondert als Download angeboten. Die Trennung wurde bewusst so vorgenommen, um dem 
experimentellen Benutzer die Möglichkeit zu geben, allein mit der Excel-Arbeitsmappe 
grundlegende Kenntnisse zu erlangen. Die wichtigsten Vorgaben und stichwortartigen 
Richtlinien sind in der Excel-Arbeitsmappe auf der Benutzeroberfläche integriert und leiten 
auch den unbelesenen Anwender zu schnellen Ergebnissen. Der interessierte Benutzer, der 
etwas tiefer in die Thematik einsteigen will, erhält die Möglichkeit die Bedienungsanleitung 
herunterzuladen und hierdurch jeden einzelnen Aspekt der Parameter nachzuvollziehen. 
Somit sind deutlich präzisere Ergebnisse möglich, die aber für den experimentellen Nutzer gar 
nicht von Bedeutung sind. Die Bedienungsanleitung kann entweder über den Link im 
nächsten Abschnitt selbst runtergeladen oder im Anhang C dieser Arbeit eingesehen werden. 


5.2 Einrichten der Webseite 


Um den Gestaltungsideen gerecht zu werden eignete sich für den Zweck der 
Webseitenerstellung die einfache Programmierung in html. Html ist eine weit geläufige 
Auszeichnungsform von Internetseiten. Es finden sich daher zahlreiche Literaturquellen, die 
die Webseitencodierung behandeln. Die Grundlagen der Programmierung, die für die erstellte 
Webseite ausreichend waren, entstammen Haunschild 2017. Da der Schwerpunkt dieser 
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Arbeit nicht auf der Programmierung einer Website liegt, sondern diese nur die Ergebnisse 
zur Verfügung stellen soll, wird auf eine detaillierte Beschreibung des Programmiervorgangs 
verzichtet. Der Quelltext der Webseite kann bei Interesse in Anhang D eingesehen werden 
und ist aufgrund des geringen Umfangs mithilfe von Haunschild 2017 zügig 
nachzuvollziehen. Die Oberfläche der Internetseite ist in Bild 5.1 dargestellt. 


Diederich Method 


Excel-tool for calculating lift distributions and other significant values for wing design purposes 


T User Manual "Diederich Method Excel-file”.pdf" Size: 72 KB" 





& Diederich Method Excel-file".xis " Size: 700 KB" 











Bild 5.1 Oberfläche der erstellten Webseite 


Die Webseite kann später als Subdomain von Professor Scholz über folgenden Link 
aufgerufen werden. 





http://Diederich.ProfScholz.de 
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6 Zusammenfassung 


Es ist eine Excel-Tabelle entstanden, die insbesondere Studierenden den leichten Einstieg in 
die Thematik der Tragflügel-Aerodynamik bietet. Neben den grundlegenden Inhalten der 
Diederich-Methode, sind noch weitere interessante Berechnungen mit in die Arbeitsmappe 
aufgenommen worden. Im Output können so auch von der Auftriebsverteilung unabhängige 
Ergebnisse für Atmosphäre oder Flügelgeometrie abgelesen werden, die für ein weites 
Themenspektrum von Interesse sein können. Der einfache Zugang über die erstellte Webseite 
ermöglicht einen stetigen Bezug zu den Daten. Die komprimierte, aber bedachte Anzahl an 
Hinweisen und die Bereitstellung der Bedienungsanleitung führen jeden Interessierten schnell 
zum gewünschten Erfolg. Die Bereitstellung der Ergebnisse in englischer Sprache ermöglicht 
es auch ausländischen Interessenten die mit der Excel-Tabelle angebotene Diederich-Methode 
zu verwenden. 
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Anhang A 


Gleichungen für die ISA 


Equations for the troposphere 


Troposphere from Om=0ft to 11000 m = 36089 ft (geopotential height) 


T=T,-L-H H geopotential height 
T, = 238.15K 
L = 0.0065 K/m = 6.5 K/km = 1.9812: 10° K/ft 
ö-2=(1-k,.H)"  p, = 101325 Pa= 1013.25hPa = 1.01325 bar 
Po 
o=2=(1-k:H)  p, = 1.225 kg/m? 
Po 
k, = 2.2558: 10° lm = 0.022558 1/km = 6.8756 ° 10% 1/f 


Equations for the Stratosphere 
Stratosphere from 11000 m = 36089 ft to 20000 m = 65617 ft (geopotential height) 
T=T, = 216.65 K = -56.5 °C = const 


o p ö p eki#-Hr) 


H, = 11000 m= 11 km = 36089 fi 

k, = 1.537688 10* 1/m = 0.157688 1/km = 4.80634 10° 1/ft 
o, = 0.297070 

Pr = 0.3639 kg/m? 

ö, = 0.223356 

Pr = 22632 Pa= 226.32 hPa = 0.22632 bar 


Bild A.1 Tropos- und Stratosphärengleichungen nach der ISA (Scholz 2012) 
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speed of sound a=,/y-R-JT MR = 20.0468 UK: 
Ss 

J 
R= 287.053 —— y=14 

kg-K 


a, = 340.294 m/s = 1225.06 km/h = 661.48 kt 


dynamic viscosity 1 ßB, > 1.458: 106 —E_ 
T+S ö m-s-/K 
k 
1. = 1.7894 10° S = 1104K 
m-s 
kinematic viscosity v= 2 v = 1.4607 10° I 
p s 
relative density 4 
Po 
. _Pp 
relative pressure ö=— 
Po 
j T 
relative temperature = T. 
0 
J 
equation of state for a perfect gas P-RT R = 287.053 —— 
p kg-K 
5, 
o 


Bild A.2 Allgemeine Gleichungen für die Atmosphärendaten nach der ISA (Scholz 2012) 
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Anhang B 


Programmcodes der in Visual Basic erstellten 
Funktionen 





(Allgemein) “| [C_two 

















Function C_two(x As Double) As Double > 
C_two = 0.000003 * (x * 5) - 0.000094 * (x * 4) + 0.0010118 * (x * 3) - 0.0015527 * (x * 2) - 0.1042267 * x + 1.001183 
End Function 
Function C_one(x As Double) As Double 
C_one = -0.0000264 * (x * 3) - 0.000564 * (x * 2) + 0.0556219 * x - 0.0054559 
End Function 
Function C_three(x As Double) As Double 
C_three = -0.0000506 * (x * 3) - 0.0007669 * (x * 2) + 0.0490423 * x - 0.002549 
End Function 
Function C_four(x As Double) As Double 
C_four = 0.0000034 * (x * 4) - 0.0000064 * (x * 3) - 0.0037242 * (x * 2) + 0.0913507 * x + 0.0008776 
| End Function 
Function F(Factor_5 As Double, x As Double) As Double 
Dim yl, y2, y3, y4, y5, y6 As Double 
Dim dyl, dy2, dy3, dy4, dy5, dy6 As Double 
Dim slopel, slope2, slope3, slope4, slope5 As Double 
Dim interpl, interp2, interp3, interp4, interp5 As Double 
= -45 
= -30 
0 
30 
45 
60 
dyl = -3.553922 * (x * 6) + 12.825226 * (x * 5) - 22.261029 * (x * 4) + 21.497862 * (x * 3) - 11.867439 * (x * 2) + 0.515369 * x + 2.143885 
dy2 = -0.131119 * (x * 4) + 1.329643 * (x * 3) - 3.354167 * (x * 2) + 0.027496 * x + 1.798392 
dy3 = -17.156863 * (x * 6) + 41.534691 * (x * 5) - 37.075792 * (x * 4) + 14.199146 * (x * 3) - 2.736741 * (x * 2) + 0.062062 * x + 1.274186 
dy4 = -11.538462 * (x * 5) + 25.670163 * (x * 4) - 22.38345 * (x * 3) + 6.815559 * (x * 2) + 0.433846 * x + 0.897133 
dy5 = 2.124183 * (x * 6) - 15.667421 * (x * 5) + 28.172448 * (x * 4) - 22.782548 * (x * 3) + 6.031885 * (x * 2) + 1.406838 * x + 0.64797 
dy6 = 2.369281 * (x * 6) - 25.53733 * (x * 5) + 50.684703 * (x * 4) - 41.851819 * (x * 3) + 13.047359 * (x * 2) + 0.677491 * x + 0.549916 
slopel = Factor_5 - yl 
slope2 = Factor_5 - y2 
slope3 = Factor_5 - y3 
slope4 = Factor_5 - y4 
slope5 = Factor_5 - y5 
interpl = (dy2 - dyl) * (slopel / 15) + dyl 
interp2 = (dy3 - dy2) * (slope2 / 30) + dy2 
interp3 = (dy4 - dy3) * (slope3 / 30) + dy3 
interp4 = (dy5S - dy4) * (slope4 / 15) + dy4 
interp5 = (dy6 - dy5) * (slope5 / 15) + dy5 
If Factor_5 < yl Then 
F = dy 
lseIf Factor_5 < y2 Then 
interpi 
Elself Factor_5 < y3 Then 
F = interp2 
Elself Factor_5 < y4 Then 
F = interp3 
ElseIf Factor_5 < y5 Then 
F = interp4 
ElseIf Factor_5 < y6 Then 
F = interp5 
Else 
F= dy6 
End If 
End Function 
I v 
|=J= < > 





Bild B.1 Ausgabefunktionen für die Faktoren C, - C, und die Pfeilungswinkelkorrekturfunktion 
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(Allgemein) v| /cLa r 








Function CLA(Factor_3 As Double, x As Double) As Double 





Dim yl, y2, y3 As Double 

Dim dyl, dy2, dy3 As Double 
Dim slopel, slope2 As Double 
Dim interpl, interp2 As Double 


yl = 1000000 
y2 = 10000000 
y3 = 100000000 


dyl = -0:.1547 * (x * 2) - 1.093 * x + 0.9 
dy2 = -0.9966 * (x * 2) - 0.6552 * x + 0.9504 
dy3 = -1.6667 * (x * 2) - 0.3822 * x + 0.9656 


slopel = Factor_3 - yl 
slope2 = Factor_3 - y2 


interpl = (dy2 - dyl) * (slopel / 9000000) + dyl 
interp2 = (dy3 - dy2) * (slope2 / 90000000) + dy2 


If Factor_3 < yl Then 
CLA = dyl 
ElseIf Factor_3 < y2 Then 
CLA = interpl 
ElseIf Factor_3 < y3 Then 
CLA = interp2 


Else 
CLA = dy6 
End If 
End Function 
v 
= < > 


Bild B.2 Ausgabefunktion für das Auftriebsgradientenverhältnis 
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Anhang C 


Online verfügbare Bedienungsanleitung der Excel- 
Arbeitsmappe 


User Guide for the Diederich Method Excel-File 


The following instructions are a compact view over the use and understanding of the 
Diederich Method Excel-file. This teaching tool can be used to calculate lift distributions for 
straight tapered swept wings. In addition several other values, useful for aircraft design 
purposes, can be determined. For in depth understanding of the calculation, read through 
the pages 467-476 in Synthesis of Subsonic Airplane Design by Egbert Torenbeek 
(available from: https://bit.Iy/3m&KlIIV). 


The Main Page 


The main page basically covers all you need for working with the Excel-file. The input fields in 
blue and the output fields in green are placed on the left side. The diagram in the middle is 
divided into two areas. The upper area shows the lift distributions over a half span. The 
curves are assigned by name in the legend, located on the right of the diagram. Below the lift 
distributions there is a visualization of the half span corresponding to the parameters stated 
in the input. On the right side there is a "How to use" field which covers the most important 
guidelines and is to be used additionally to this user guide. Plot tables are available at the 
bottom ofthe main and any other page, where they can be visited if necessary. 


The Input 


The input is divided in two sections. The General Input is the input with the most affecting 
parameters on the displayed diagrams. The restrictions for the general input are stated in 
the “How to use” guide on the right and therefore not mentioned furthermore. The 
Additional Input covers the airfoil and atmosphere data and gives an opportunity to 
differentiate between different airfoils for the tip and for the root chord. If the airfoil data is 
not given, see Theory of Wing Sections by Abbott and Doenhoff from 1959. The definition of 
the trailing edge angle of the airfoil is given in Figure 1 on the next page. The airfoil for the 
default setup of the additional input is taken from NACA 63;-615, also listed in Theory of 
Wing Sections by Abbott and Doenhoff. For determining the correct airfoil data, the 
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Reynolds number is given in the output section for the chosen flight altitude and the chosen 
mean chord length. Note: To perform the calculations, it is not necessary to know the type 
of airfoil. It is sufficient to know (or to guess) the numbers. 


The Output 


The Output is divided into three sections. The Comparison Output shows the deviations 
between the calculated lift distribution for the chosen wing and the comparative 
distributions. The elliptic distribution stands for minimum induced drag while the triangular 
distribution provides a relatively small root bending moment. The Mason-Curve is something 
in between and combines both benefits to an optimum distribution. The distribution with 
the smallest number stated has the most similarity to your calculated distribution and 
therefore indicates the given characteristics. The Additional Output provides the 
intermediate and additional results that are potentially of interest to the user. The 
Calculated Values Output shows numerous intermediate results of the calculation of the lift 
distributions in a bundled manner. Results can be tracked here, and, in the case of manual 
calculation, these are available to the user. Most of these results can also be found in the 
additional worksheets with the respective product graphs. 


The Additional Worksheets 


The additional worksheets are only for the more interested users and include the individual 
sub-areas of the calculation of the lift distributions and therefore give brief insights into the 
topic. For a full understanding, the literature mentioned in the beginning should be noted. 





an — 


Figure 1 Definition of the trailing edge angle of the airfoil pr; 
(Hoak, D.E.: USAF Stability and Control Datcom) 


